Nonuniform current and spin accumulation in a 1 \im thick n-GaAs channel 
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The spin accumulation in an n-GaAs channel produced by spin extraction into a (Ga,Mn)As contact is 
measured by cross-sectional imaging of the spin polarization in GaAs. The spin polarization is observed in a 
1 /im thick n-GaAs channel with the maximum polarization near the contact edge opposite to the maximum 
current density. The one-dimensional model of electron drift and spin diffusion frequently used cannot explain 
this observation. It also leads to incorrect spin lifetimes from Hanle curves with a strong bias and distance 
dependence. Numerical simulations based on a two-dimensional drift-diffusion model, however, reproduce the 
observed spin distribution quite well and lead to realistic spin lifetimes. 



Spin injection from a ferromagnetic contact into a 
semiconductor is a fundamental prerequisite for many 
spintronic device J^^. Moreover, a reasonably large spin 
lifetime in the semiconductor is necessary for any useful 
application. The value of the spin lifetime is often ex- 
tracted from Hanle-curves obtained from non-local volt- 
age measurements or by optical mean^^"^''^''^ using a one- 
dimensional spin drift-diffusion model. Parasitic contri- 
butions such as dynamic nuclear polarization at low tem- 
peratures have to be avoided by fast periodic magnetiza- 
tion reversal of the injecting contactP. Here we show that, 
in addition, the electric field near the injecting contact 
strongly influences the spin density distribution around 
it. As a result a one-dimensional picture of drift and 
diffusion along the channel leads to incorrect spin life- 
time values. This becomes obvious by comparing fitting 
results of Hanle data from a one- and a two-dimensional 
model. For our studies we employ a cross-sectional imag- 
ing method that allows to probe the two-dimensional 
spin polarization distribution even below the contacts^iSj 
and compare the experimental results to two-dimensional 
drift-diffusion simulations. 

The measurement principle and sample geometry are 
sketched in Fig. [l] The layers are grown by molecu- 
lar beam epitaxy on a semi-insulating GaAs(OOl) sub- 
strate. The 1 |J.m thick n-GaAs channel for the electron 
transport and spin accumulation has a doping density 
of 2.7 • 10^® cm^'^. Spin injecting contacts are formed 
by a 50 nm thick layer of the diluted magnetic semi- 
conductor GagsMusAs, a 8 nm thick layer of n++-GaAs 
doped up to 5 • 10^® cm~'^ and a 15 nm thick n++ n 
transition layer. Together these layers form a spin Esaki 
diode that enables tunneling of spin-polarized electrons 
between (Ga,Mn)As and the GaAs channel. The de- 
tailed layer stack and the contact patterning are de- 
scribed elsewher^lfil, xhe fabricated lateral device has 
a 50 |J.m wide mesa channel, one 4 [a.mx50 |J.m inject- 
ing contact and two large reference contacts located at a 
distance of 300 |J,m to the Ga(Mn,As) contact. Finally, 
the sample is cleaved along the [110] direction across the 
mesa channel and the ferromagnetic pad, reducing the 
contact area to 4 |xmx40 |J.m. The cleaving process fur- 
thermore exposes the GaAs(l-lO) surface and enables a 
direct optical access to the n-GaAs channel (see Fig. [ij. 
For the optical measurements the sample is mounted in a 
He flow cryostat. The cryostat itself is mounted on top of 




FIG. 1. Geometry of the sample and measurement principle. 



a nano positioner enabling two-dimensional scans in the 
xy-plane. The z component of the electron spin polariza- 
tion (that is, the component along [1-10], which is a mag- 
netic easy axis of the Ga(Mn,As) contact) in the n-GaAs 
channel is detected via the polar magneto-optical Kerr ef- 
fect (pMOKE) with a spot size of 1 |im. The photon en- 
ergy of the linearly polarized laser beam was chosen close 
to the band-gap of GaAs (A = 816 nm at 9 K), where 
the specific Kerr rotation shows a maximum (not shown) . 
A square- wave bias voltage alternating between zero and 
Vb is applied between the contacts and the Kerr rotation 
is detected synchronously with balanced photo-receivers 
and a lock-in technique. This ensures that the (quasi- 
static) magnetization of the ferromagnetic contacts does 
not contribute to the Kerr signaF''^. 

Fig.[2]shows one-dimensional scans of the Kerr rotation 
along the n-GaAs channel for different applied positive 
voltages between the (Ga,Mn)As contact in the middle 
and the reference contact at the right hand side. As 
illustrated in Fig. [2] this denotes that unpolarized elec- 
trons in the channel are accelerated to the left and the 
spin accumulation is generated by an electron flow into 
the (Ga,Mn)As contact. Thus, spins are extracted from 
the n-GaAs layer. To eliminate any electro-optic back- 
ground, the measurements are performed in remanence 
after saturation along [1-10] and [-110], respectively and 
the difference between both remanent values is used as a 
measure of the spin polarization in the GaAs layer. 

The observed spin density distribution in Fig. [2] shows 
an unexpected maximum spin polarization that is located 
near the contact edge opposite to the maximum charge 
current. This peak shifts even further to the left with 
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FIG. 2. Kerr rotation along the n-GaAs channel for various 
applied voltages, representing the local spin polarization. 



increasing bias voltage. On the left hand side of the con- 
tact the spin density distribution should be solely driven 
by diffusion. However, the measured decay of the spin 
polarization on the diffusion side cannot be described by 
a single exponential decay as expected (not shown) . De- 
pending on the bias voltage, a faster decay close to the 
contact edge is observed. Since the decay constant de- 
fines the spin diffusion length Lg = V Ds • , either the 
spin diffusivity or the spin lifetime must be reduced in 
this area. A small bias dependence is still visible when fit- 
ting the curves from y-position —30 |im to —7 |im, where 
a spin diffusion length of 8.2 |j.m was extracted for the 
largest bias (9V / 750 |J.A) and 9.2 \xm for the smallest 
bias (IV / 25 M-A). On the right hand side of the con- 
tact, the drift side, the spin distribution can be described 
by a superposition of drift and diffusion. This explains 
the shortened effective decay length with increasing bias 
voltage, but it does not explain the decay of the spin 
polarization that already shows up beneath the contact 
area. As a verification of this spin density distribution, 
we also perform non-local voltage measurements on a 
similar sample with additional (Ga,Mn)As contacts on 
the drift and diffusion side. The result from the electri- 
cal detection experiment shows a good agreement with 
the optical data (not shown). 

For a better understanding of the observed spin den- 
sity distribution and its bias dependence we also perform 
Hanle-measurements at y-position —5 |J.m and —10 [a.m 
for different applied voltages, illustrated in Fig. [31 The 
curves are fitted using a one-dimensional functioiP'. As- 
suming pure spin diffusion from the contact area to the 
measurement position, the only free parameter left is the 
spin lifetime, as the spin diffusion length is well known 
from the previously shown results (see Fig. [2]). How- 
ever, this approach leads to spin lifetimes depending sig- 
nificantly on the applied bias voltage and the injection- 
detection distance (see Fig. |3| , with values changing from 
14 ns for a large voltage and small distance to 46 ns for a 
small voltage and large distance. The bias dependence is 
less pronounced when measuring further away from the 
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FIG. 3. Hanle-curves for different applied bias voltages, mea- 
sured at y-position —5 \im (top) and —10 \im (bottom). Fits 
using a one-dimensional model yield different values of the 
spin lifetime which are shown in TABLE I. 



contact, but still a factor of two difference is observed. 
However, even the IV-data with a current of 25 |J.A yield 
different values at both positions (42 and 46 ns). 

To shed light on this behavior we perform numerical 
simulations of the two-dimensional electron drift and spin 
diffusion. The simulation is done in two steps: first the 
current density distribution is simulated in the n-GaAs 
channel that in principle depends on the interface resis- 
tance of the Esaki diode and the channel resistance. Both 
values are determined by electrical measurements, where 
the voltage drop across the interface is detected via the 
second reference contact. At 9 K, the channel's conduc- 
tivity is 1500 S/m, resulting in an electron mobility fi 
of 0.35 m^/(s • V) in the n-GaAs channel. In the simula- 
tion, the Esaki diode is represented by a 10 nm thick layer 
with a conductivity of 0.025 S/m, matching with the ex- 
perimentally determined resistance for this bias voltage. 
Based on the result for the current density distribution, 
the spin density distribution is calculated by assuming 
a spin current across the interface that is proportional 
to the current density. The generated spin accumulation 
then spreads due to diffusion and the previously calcu- 
lated drift in this sample geometry until a steady state 
is found. A derivation of the used differential equations 
and boundary conditions can be found in Ref. 2. The 
spin diffusivity in the simulation is defined by the spin 
diffusion length, extracted from the measurements shown 
in Fig.[2]and the spin lifetime (39ns), estimated by Hanle- 
measurements (see below). 

Fig. [4] illustrates the simulated spin density distri- 
bution for a large positive applied bias voltage at the 
(Ga,Mn)As contact, which reproduces the observed spin 
density distribution quite well: the pronounced nonuni- 
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FIG. 4. Two-dimensional drift-diffusion simulation of the spin 
accumulation in a 1 \im thick n-GaAs channel. Black lines 
represent the local current direction. 

form electron drift beneath the contact area shifts the 
polarization peak towards the diffusion side. This shift 
increases with applied bias voltage similar to our exper- 
imental results (see Fig. [2]) . The simulation also directly 
explains the extremely fast decay of the spin polarization 
on the right hand side of the contact, where the decay 
already begins beneath the contact area. 

An important point to note is the nonuniform current 
density on the diffusion side and underneath the contact 
area. As illustrated in Fig. |4] by the black lines repre- 
senting the local current direction, the n-GaAs channel 
on the diffusion side can be divided into a bottom chan- 
nel where electrons are flowing to the left away from the 
contact and a top channel with a back-flow towards the 
contact. This nonuniform current density decays rapidly 
with increasing distance but cannot be neglected within 
the first micrometer on the diffusion side. Due to this 
inhomogeneous current density it is not sufficient to de- 
scribe the electron spins solely by diffusion on this side 
of the channel, which was done before when fitting the 
Hanle-curves. fn particular, the electric field beneath 
the contact area generates a positive drift towards the 
measurement position. This additional drift widens the 
Hanle-curves and therefore the fit provides lower spin life- 
times. 

This effect becomes apparent when the Hanlc data is 
fitted with our two-dimensional model where the bias 
dependent spin density distribution and the nonuniform 
drift is taken into accountf^i': 
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TABLE I. Spin lifetime from numerical simulations based on 
ID and 2D drift-diffusion models at two different distances, 
y, from the contact edge and different bias voltages'^. 

A spin lifetime of 38 ns is extracted for the 9 V data in 
contrast to 14 ns for the one-dimensional model. Thus, 
the strong bias dependence of the spin lifetime from the 
one-dimensional model can be primarily attributed to the 
neglect of the nonuniform current. However, a variation 
of the spin lifetime with bias voltage and contact distance 



of about 20% remains. 

The dominant spin relaxation mechanism for n- 
doped GaAs above the metal-insulator transition is the 
Dyakonov-Perel mechanisirP^]. This mechanism strongly 
depends on the electron energy in the conduction band, 
for instance on the doping density or the temperature 
(in case of charged impurity scattering oc 7^-3/2^^ 
According to this formula, a temperature increase up to 
2 K would be sufficient to explain the observed behavior. 
Besides the temperature, also the electric field has an in- 
fluence on the spin lifetime, as has already been shown 
for optically pumped spins in 200( P^P^ . The reason is 
the increased energy of the electrons when accelerated in 
electric fields leading to an increased momentum relax- 
ation time Tp and thus to a more efficient spin relaxation 
due to the Dyakonov-Perel mechanism^. A detailed cal- 
culation of this effect was reported by Beck et 

From the one-dimensional scans and the simulation it 
is evident that electric fields are present around the con- 
tact area. Even for the IV data, the electric field in the 
n-GaAs channel is in the order of 10 V/cm, that, accord- 
ing to Refs. 12 and 13 should already reduce the spin life- 
time. Since strong electric fields are only located beneath 
the contact area and the electric field intensity decreases 
rapidly on the diffusion side, larger spin lifetimes and a 
less pronounced bias dependence are observed when mea- 
suring further away from the contact. 

In summary, it is shown that a nonuniform current 
density in the n-GaAs channel affects the spin density 
distribution. The pronounced electric field beneath the 
contact area generates a drift towards the diffusion side 
and therefore widens Hanle-curves when measured on 
this side of the contact. As a consequence, extracting 
the spin lifetime with a one-dimensional model assuming 
pure spin diffusion leads to a strong bias and contact dis- 
tance dependence. Hence a two-dimensional model has 
to be used to yield correct spin lifetimes. The remaining 
variations with distance and bias may be due to Joule 
heating or the presence of electric fields around the con- 
tact area. As a final note, the nonuniform current density 
might also be the origin of offset-signals in non-local volt- 
age measurements and directly shows that the non-local 
voltage is influenced by resistance changes of the injector 
contact. 
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